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ABSTRACT 

We present results showing that the multiple system DENIS-P J020529. 0-115925 is likely to 
be a triple system of brown dwarfs. The secondary of this previously known binary system shows 
a clear elongation on six images obtained at six different epochs. Significant residuals remain 
after PSF subtraction on these images, characteristic of multiplicity, and indicating that the 
secondary is probably a double itself. Dual-PSF fitting shows that the shape of the secondary 
is consistent with that of a binary system. These measurements show that the probability that 
DENIS-P J020529. 0-115925 is a triple system is very high. The photometric and spectroscopic 
properties of the system are consistent with the presence of three components with spectral types 
L5, L8 and TO. 

Subject headings: Stars: individual (DENIS-P J020529. 0-115925) - Stars: low mass, brown dwarfs - 
Binaries: visual, triple system - Techniques: high angular resolution 
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1. Introduction 

Multiple systems are important tests for the 
models of formation and evolution of very low 
mass stars and brown dwarfs. The binary fraction 
reported recently by Bouy et al. (2003); Burgasser 
ct al. (2003); Close et al. (2003); Gizis et al. (2003); 
Martin et al. (2003) for very low mass stars and 
brown dwarfs (between 10^15%) cannot be well 
reproduced by the so-called "photo-evaporation" 
and "embryo-ejection models", which predict 
a much lower binary frequency (<5%; see e.g 
Kroupa & Bouvier 2003, for a discussion on the 
different scenarii of formation mentioned here). 
On the other hand the "star-like" model cannot 
explain the observed distribution of separation (all 
less than 20 A.U, with a peak around 4~8 A.U), 
and of mass ratios (with a strong lack of multi- 
ple systems with large mass ratios). The overall 
process of formation of brown dwarfs is still not 
well understood. The existence itself of an old 
triple system made of brown dwarfs would give 
extremely important and new constrains on these 
models, since no such object has been observed 
yet. 

In section 2, we will present a summary of the 
known properties of DENIS-P J020529.0-1 15925. 
In section 3 we will describe the observations. In 
section 4, we will present the analysis of the data 
leading to the discovery of a possible third com- 
ponent, and in section 5 we will discuss about the 
triple system and its properties. 

2. DENIS-P J020529.0-115925 

DENIS-P J020529.0-1 15925 is a L5 field very 
low mass dwarf (Martin et al. 1999). From 
their photometry and spectroscopy, Geballe et al. 
(2002) classify the object MS cX L5.5±2 dwarf. It 
is classified as L7 in the Kirkpatrick et al. (2000) 
scheme. We adopt a spectral class of L5. It was 
discovered by Delfossc et al. (1997) and first re- 
solved as a binary by Kocrncr et al. (1999) using 
Keck images. This object has been well studied 
and observed, and is reported in several surveys 
(DENIS, and 2MASS, as 2MASSW J0205293- 
115930). Several authors (Delfosse et al. 1997; 
McLean et al. 2001; Burgasser et al. 2003) report 
methane absorption in their spectra, which im- 
plies a mass below the sub-stellar limit and an ef- 
fective temperature less than 1800 K, as stated by 



Schweitzer et al. (2002). Tokunaga & Kobayashi 
(1999) report similar absorption features in their 
spectra but attribute it to H 2 rather than CH4. 
Martin et al. (1997) reported a non-detection of 
lithium absorption from high resolution optical 
spectra, and inferred a lower limit on the mass 
of DENIS-P J020529.0-115925 A of 60 Jupiter 
masses. Its distance (19.76±0.57 pc) and proper 
motion (437.8±0.8 mas/yr with PA=82.8±0.1°) 
have been measured via trigonometric parallax 
(Dahn et al. 2002). Basri et al. (2000) also mea- 
sured its rotational velocity (22±5 km/s), and 
from the Cs I and Rb I absorption, they estimated 
its temperature to be T e// = 1700~1800 K. 

3. Observations 

We observed DENIS-P J020529.0-115925 using 
the Hubble Space Telescope (HST) at six differ- 
ent epochs. The observations occurred between 
October 2000 and December 2003 during HST 
Cycles 8, 9, 10, 11 and 12 (program GO8720, 
P.I. Brandner, and programs G09157, G09345 
and G09968, P.I. Martin). DENIS-P J020529.0- 
115925 was observed with the Planetary Cam- 
era of WFPC2 (Biretta 2002), in the F675W and 
F814W filters (GO8720) and F606W and F814W 
filters (G09157, G09345 and G09968). Table 1 
gives a log of all the observations we used for this 
study. The target is very red, and the observations 
in the F814W filter were more sensitive than the 
one in the F675W or F606W filters, despite the 
shorter exposure times and the lower quantum ef- 
ficiency of WFPC2 at longer wavelengths. The 
possible third component does not appear clearly 
on the F675W and F606W images. For these rea- 
sons, the following analysis is made essentially in 
the F814W filter. The high angular resolution of 
the WFPC2-PC (tf'0455 pixel scale) allowed us to 
resolve easily the 0'.'3 secondary, and to discover 
that it is elongated. 

4. Data analysis 

We processed the data in two steps: a point 
spread function (PSF) subtraction on each compo- 
nent of the binary to look for residuals indicating 
the presence of a third companion, and a dual-PSF 
fitting of the secondary for confirmation. 
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4.1. PSF Subtraction and Residuals 

Taking advantage of the extremely stable PSF 
of HST/WFPC2, we performed a PSF subtraction 
on the six images of DENIS-P J020529.0-115925 
using 9 different unresolved objects from programs 
GO8720 (P.I. Brandner, Bouy et al. 2003) and 
G08581 (P.I. Reid, Gizis ct al. 2003), and 1 syn- 
thetic PSF obtained with Tiny Tim (Krist & Hook 
2003) as reference PSF stars. Figure 1 shows that 
on 4 of the 6 images, the primary is extremely well 
subtracted, while strong residuals appears clearly 
after subtraction of the secondary, indicating the 
presence of a third component. The results with 
the 9 other reference PSF stars are similar to that 
presented in Figure 1. 

Using this technique, we tried to define a cri- 
terion that would allow to determine quantita- 
tively if the secondary is a multiple system or 
not. We performed a detailed PSF analysis of 
all the WFPC2 images of programs G08581 (P.I. 
Reid Gizis et al. 2003). This large sample has the 
advantage to be homogeneous (all targets were ob- 
served with the F814W filter of the WFPC2/PC), 
and all targets are very low mass stars and brown 
dwarfs with spectral properties (and therefore 
PSF) similar to DENIS-P J020529.0-115925. 

The relative intensity of the residuals, defined 
as the integrated intensity of the residuals after 
PSF subtraction, divided by the integrated inten- 
sity of the object as expressed in equation 1, ap- 
pears to be a powerful criterion to identify binary 
candidates. 

n.i=^ (i) 

i=i 

where: n is the number of pixels, J-r{i) is the flux 
in pixel i after PSF subtraction, and To{i) is the 
flux in pixel i before PSF subtraction. 

Figure 2 shows that 1Z.1 is very low and 
very stable for unresolved objects, while it is al- 
ways higher than ^3-cr above the median value 
in the case of multiple systems. This tech- 
nique found easily and automatically all the 
multiple systems present in the sample (see 
Bouy et al. 2003; Gizis et al. 2003), except 
one (2MASSW J0856479+223518). This lat- 
ter multiple system is indeed a good example 



to illustrate the limitations of this technique. 
2MASSW J0856479+223518 is a very close binary 
(5=0>! 1) with a relatively large difference of mag- 
nitude (AMag=2.8 mag, see Bouy et al. 2003). 
The relative intensity of the residuals after PSF 
subtraction of the primary, corresponding roughly 
to the ratio of the intensity of the secondary over 
that of the primary, is therefore of the order of a 
few percent only (AMag=2.8 mag is equivalent to 
a flux ratio of ~0.07). The relative inten- 

sity of the residuals after PSF subtraction is there- 
fore a good method to find multiple systems with 
moderate differences of magnitude. Fortunately, 
it is the case of the possible third component of 
DENIS-P J020529. 0-1 15925, as shown below. 

Using this property, we compare the relative in- 
tensity of the residuals after PSF subtraction on 
the secondary of DENIS-P J020529.0-115925, us- 
ing the same library of PSF stars. Figure 2 shows 
the results. The residuals at six different epochs 
are all higher than that of unresolved objects and 
unresolved companions of known binaries of the 
same sample, but very similar to that of other 
multiple systems. This figure shows clearly that 
DENIS-P J020529.0-115925B is very likely to be 
a binary system itself. Figure 2 shows also that 
even for the two images where the residuals do 
not appear clearly by eyes (see Figure 1), the 1Z.T 
is slightly higher than in any unresolved object. 

One should note that 1Z.I strongly depends on 
the signal-to-noise ratio (S/N). At low S/N, the 
value of must indeed increase due to the in- 
creased noise in the target. Figure 2 shows that all 
our images were obtained at sufficiently high S/N 
(always greater than 100), so that 1Z.I provides a 
robust statistical test of multiplicity. 

4.2. PSF fitting 

In order to confirm that these residuals are con- 
sistent with the presence of a third component, we 
performed a dual-PSF fitting, on the secondary 
only, using a custom-made program described in 
Bouy et al. (2003). Briefly, and as presented in 
Figure 3, the PSF fitting routine builds a model 
binary using each of the ten different PSF stars 
mentioned above, and then performs a non-linear 
PSF fit of the observed image to determine the 
best-fit values for the three free parameters: sepa- 
ration, position angle and flux ratio. The average 
of the ten results gives the final values for the pa- 
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rameters. Figure 3 shows that the residuals after 
this dual-fit are much better than after the singlc- 
PSF subtraction described in the previous section 
and in Figure 1. Moreover, the stability of the re- 
sults for the ten different PSF stars at each of the 
six epochs shows that the results are reliable and 
robust, and that the third component hypothesis 
is consistent with the shape of the elongated PSF 
and highly probable. 

Table 2 gives the flux ratios between the three 
components of the multiple system estimated with 
this method, as well as the corresponding differ- 
ences of magnitude. According to the most recent 
DUSTY models of Chabrier et al. (2000), these 
differences of magnitudes indicate that the three 
components must have similar masses. A differ- 
ence of magnitude of 1.5 mag in the F814W cor- 
responds indeed to a mass ratio of 75, 85 and 95% 
at respectively 0.5, 1.0, and 5 Gyr. 

5. Discussion 

5.1. Properties of the triple system 

The present analysis of the high angular reso- 
lution images indicates that DEN1S-P J020529.0- 
115925 is very likely to be a triple system. 

Figure 5 shows the Mi c vs Spectral Type re- 
lation for all the objects reported by Dahn et al. 
(2002) (see their Tables 1, 2 and 3). Assuming 
differences of magnitude in Ic equal to that re- 
ported in Table 2 for the F814W filter (F814W 
filter is close to I c , Biretta 2002), and the DENIS 
I magnitude 1 of the unresolved objects, one can 
estimate the spectral types of the three compo- 
nents. DENIS-P J020529.0-1 15925 A is consistent 
with a L5.5 dwarf, in good agreement with the 
measurement reported by Martin et al. (1999) and 
Geballe et al. (2002), and showing that the pri- 
mary would be dominating the optical spectrum. 
B and C would be consistent with — L8 and —TO 
dwarfs respectively, in the Geballe et al. (2002) 
classification scheme. 

Several authors report the detection of methane 
absorption in the infrared spectrum of DENIS- 
P J020529. 0-1 15925 (Delfosse et al. 1997; McLean 
et al. 2001; Burgasser et al. 2003). Tokunaga & 
Kobayashi (1999) also observed this absorption 
feature in the spectrum but attribute it to H 2 

1 Idenis is very close to the Icousins (Delfosse 1997) 



rather than methane. Burgasser et al. (2002) note 
that this feature is weak and variable, and con- 
sider that it does not constitute a clear detection 
of methane. We note that if real, this feature could 
be related to the presence of L8 and TO compan- 
ions, and its variability to some weather effects, as 
already observed in other late L and T dwarfs by 
Enoch et al. (2003). Dahn et al. (2002) showed 
that the absolute J and K-band magnitudes of 
early T-dwarfs are similar to that of late-L dwarfs, 
so that the contribution of B and C in the near- 
infrared can be larger than in the optical. 

According to the DUSTY models (Chabrier 
et al. 2000), and assuming an age between 1 
and 10 Gyr, Figure 6 shows that the absolute 
Mi c magnitude of DENIS-P J020529.0-115925 
A corresponds to an effective temperature be- 
tween 1700-1900 K, while B and C range be- 
tween 1 550-1 800 K. These temperatures are con- 
sistent with the value reported by Basri et al. 
(2000) for the unresolved system (1700-1800 K), 
and show that all components appear to be 
clearly substcllar. According to the DUSTY mod- 
els, the stellar /substcllar limit is indeed around 
2000 K at 10 Gyr and 2180 K at 1 Gyr, therefore 
warmer than any of the components of DENIS- 
P J020529. 0-115925. 

Finally, the proper motion of the object (—438 
mas yr _1 , Dahn et al. 2002) and the presence 
of these strong residuals at six different epochs 
spread over three years allow us to rule out defini- 
tively the eventuality of a coincidence with some 
background object. 

5.2. Dynamical Properties 

The separation between the primary and the 
secondary changed from — 0'.'390 to 0'.'270 between 
October 2000 and December 2003, while the sep- 
aration between the second and the third compo- 
nent is contained between — 0'.'075< Sbc <0'/055 
(see Table 3). As stated by Harrington (1968) and 
then Szcbchcly & Zare (1977) in their analytical 
study, triple systems with moderate eccentricity 
and equal mass components are stable for ratios 
between the semi- major axes of the outer (a 2 ) 
and the inner orbits (a\) greater than a 2 /ai >3.2. 
Assuming that the orbits of DENIS-P J020529.0- 
115925 components have moderate eccentricities, 
and considering that the 3 components must 
have similar masses as explained above, DENIS- 
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P J020529.0-1 15925 would fit above the stability 
criterion, with a ratio between 3.6< 02/01 <7.1. 
The presence of a third component at the positions 
reported here is therefore dynamically possible. 

The estimate of the separation corresponds to a 
semi-major axis of ~1.2 A.U (~0'.'075 at 19.7 pc). 
Corrected for a statistical factor of 1.26 as ex- 
plained in Fischer & Marcy (1992), it leads to a 
"best-guess" semi- major axis of 1.9 A.U. Accord- 
ing to Kepler's Third Law (Kepler 1609) and as- 
suming a total mass of ~0.1 M Q , the correspond- 
ing period, is ~8 years. 

The dual-PSF fitting indicates that over the six 
epochs the separation between B and C is bounded 
between 0'.'053-0'.'074, while the position angle is 
bounded between 63° and 109° (see Table 3). The 
raw results of the PSF-fitting give a motion of C 
around B that cannot be orbital in nature (P. A 
does not vary monotically, as shown in panel a) 
of Figure 4), but at such very short separations, 
the possible confusion between the primary and 
the secondary during the PSF-fitting yields to an 
ambiguity of ±180° in the measurements of the 
position angle. Figure 4 illustrates two cases for 
which two epochs have been switched by 180°, giv- 
ing possibly Keplerian orbital motions. The corre- 
sponding period (about half of the orbit in about 
3 years) is in good agreement with that reported 
above (8 years). 

With large and hardly assessable uncertainties 
and with 180° ambiguity on the P.A, we consider 
that the dual-PSF fitting is essentially a valuable 
sanity check showing that the elongation of the 
secondary is consistent with a binary. It gives 
a very rough and tentative estimate of the pho- 
tometric and astrometric properties of the possi- 
ble third component. New observations at higher 
angular resolution, using for example HST/ACS, 
ground based adaptive optics or speckle imaging, 
should be performed in order to confirm the mul- 
tiplicity, and if confirmed to obtain more precise 
measurements of the motion of C around B. 

6. Conclusions 

We present here results of high angular reso- 
lution observations with HST/WFPC2 that allow 
us to conclude that DENIS-P J020529.0-115925 is 
very likely to be a triple system of brown dwarfs. 
PSF subtraction on the secondary at six different 



epochs show unusual residuals in comparison with 
unresolved objects. Dual-PSF fitting shows that 
the shape of the secondary can be consistent with 
that of a binary, and that the motion of the possi- 
ble third component could be Keplerian in nature. 
The configuration is consistent with a dynamically 
stable multiple system. 

Observations at higher angular resolution, us- 
ing for example the HST/ACS, ground based 
adaptive optics or speckle imaging, should allow to 
confirm if DENIS-P J020529.0-115925 is a triple 
system or not. 

This work is based on observations collected 
with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. These observations 
are associated with programs GO8720, G09157, 
G09345 and G09968. This publication makes use 
of data from the DEep Near Infrared Survey. The 
DENIS project has been partly funded by the SCI- 
ENCE and the HCM plans of the European Com- 
mission under grants CT920791 and CT940627. It 
is supported by INSU, MEN and CNRS in France, 
by the State of Baden- Wrttemberg in Germany, 
by DGICYT in Spain, by CNR in Italy, by FFwF- 
BWF in Austria, by FAPESP in Brazil, by OTKA 
grants F-4239 and F-013990 in Hungary, and by 
the ESO C&EE grant A-04-046. This publication 
makes use of data products from the Two Micron 
All Sky Survey, which is a joint project of the 
University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center /California Institute 
of Technology, funded by the National Aeronau- 
tics and Space Administration and the National 
Science Foundation. This work made use of data 
from the The Guide Star Catalogue II, which is a 
joint project of the Space Telescope Science Insti- 
tute and the Osservatorio Astronomico di Torino. 

Facilities: HST(WFPC2). 
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Table 1 
Observation log 



Filter 


Exp. Time 
[s] 


Date Obs. [UT] 

DD / MM /YYYY 


Program 


F814W 


600 


28/10/2000 


GO8720 


F675W 


300 


28/10/2000 


GO8720 


F814W 


1700 


08/07/2001 


G09157 


F814W 


1800 


21/01/2002 


G09157 


F814W 


400 


14/07/2002 


G09345 


F606W 


1600 


14/07/2002 


G09345 


F814W 


400 


10/01/2003 


G09345 


F814W 


800 


01/12/2003 


G09968 


F606W 


1000 


01/12/2003 


G09968 
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Table 2 

Flux ratios and difference of magnitudes in the F814W filter 



Components Flux Ratio AMag 



B/A 0.37 1.08 

C/A 0.25 1.50 

C/B 0.66 0.45 



Note. — Abs. Mag(A)=17.3 mag 
The measurements reported here have 
large uncertainties, and should be con- 
sidered with caution. The relative er- 
ror from one image to the other is about 
0.15 mag, but the real uncertainties can 
be higher. Measurements obtained on 
the GO8720 image. 
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Table 3 

Relative Astrometry of the three components 



Date Obs. Separation P.A 



Component B with respect to A 



01/12/2003 


0'.'287±0'.'005 


246°±1° 


01/10/2003 


0"306±0'/005 


249°±1° 


14/07/2002 


0'.'330±0'.'005 


251°±1° 


21/01/2002 


0"350±0'/005 


253°±1° 


08/07/2001 


0'.'370±0'.'005 


256°±1° 


28/10/2000 


0'.'398±0'.'005 


259°±1° 


Component C with respect to B 


28/10/2000 


0'.'073±0'.'005 


86°±3° 


08/07/2001 


0'.'050±0'.'010 


63°±5° 


21/01/2002 1 






14/07/2002 


0'.'058±0'.'010 


109°±5° 


10/01/2003 


0"063±0"008 


94°±5° 



01/12/2003 1 



Note. — Component C is barely re- 
solved. The measurements reported 
for component C and their uncertainties 
are close to the limit of resolution of 
HST/WFPC2, and should be considered 
with caution. 

x The secondary is elongated, and 
the residuals after PSF subtraction arc 
strong, but the PSF fitting does not give 
good enough results on the astrometry 
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Fig. 1. — Results of the PSF subtraction at different epochs. This figure shows surface plots of the results of 
the PSF subtraction with one of the 10 reference PSF used. The primary is well subtracted whereas much 
stronger residuals remain for the secondary (clearly noticeable in 4 of the 6 epochs, indicated by arrows), 
indicating the presence of a third component. Similar results are obtained with the 9 other reference PSF 
stars. X and Y axis are in pixels, Z axis is in counts. 
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Fig. 2. — Residuals after PSF subtraction on HST/WFPC2 images (F814W filter) of very low mass stars 
and brown dwarfs. This plot shows TZ.l as a function of the signal-to- noise ratio for the targets of program 
G08581. The crosses represent the multiple systems. The plus represent the values obtained at 6 epochs for 
DENIS-P J020529.0-115925B. A filled square and a disk represent the results for two unresolved companions 
of known multiple systems, respectively 2MASSW J1146+2000B (program G08146, P.I. Rcid) and DENIS- 
P J1441-0945B (program GO8720, P.I. Brandner). The 3-er above the median and the median values for 
the unresolved objects are represented. This figure shows that the intensity of the residuals after PSF 
subtraction appears to be a powerful indicator of multiplicity. DENIS-P J020529.0-115925B appears clearly 
different from the unresolved objects and from the unresolved companions of known multiple systems, while 
its residuals are comparable to that of known close multiple systems. 
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Fig. 3.— Results of the dual-PSF fitting on DENIS-P J020529.0-f f 5925B for the GO8720 image. As shown 
with arrows, DENIS-P J020529.0-115925B has a FWHM much wider than the reference PSF star. The 
companion appears clearly after subtraction of the model primary. The residuals are small (lower panel, on 
the left), indicating that the quality of the fit is good. They are especially smaller than that obtained with 
a single PSF subtraction (lower panel, on the right). Similar results have been obtained with the 9 other 
reference PSF stars, showing that the result of the fit is robust and reliable. X and Y axis are in pixels, Z 
axis is in counts. 12 
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Fig. 4. — Panel a), b) and c): Relative motion of C (diamonds) around B (triangle). Panel d): Relative 
motion of B+C (diamonds) around A (triangle). 1-a uncertainties are indicated, as well as the first epoch. 
Epochs are linked by a dot line in chronological order. Due to the very short separation, there is an ambiguity 
of 180° in the position angle of C around B. Panels a), b) and c) illustrate three configurations among the 
possible ones, a) raw results of the PSF fitting; b) epochs 1 and 2 from the raw results but 180° is added 
to epochs 3 and 4; c) epochs 3 and 4 from the raw results but 180° is added to epochs 1 and 2. While case 
a) cannot represent any Keplerian orbit, cases b) and c) could, with orbital periods comparable with the 
preliminary estimate we report here. Values from Table 3. 
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Fig. 5. — Mi c vs Spectral Type relation. All measurements from Dahn et al. (2002) (See their Figure 2, and 
Tables 1, 2 and 3) except for DENIS-P J020529.0-115925A, B and C, evaluated as explained in the text. The 
absolute magnitudes of DENIS-P J020529.0-115925A, B and C have been estimated using the differences of 
magnitude reported in Table 2, and the DENIS I magnitude of the unresolved system. 
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Fig. 6. — M\ c vs Effective temperature. This figure shows the Mi c as a function of the effective tem- 
perature, as given by the DUSTY models. Two isochrones (1 and 10 Gyr) are shown, together with the 
corresponding limits for sustained hydrogen burning. The position of the three components is represented by 
grey boxes, assuming ±100 K uncertainty on the effective temperatures of the three components of DENIS- 
P J020529. 0-115925, and including the uncertainty on the distance for the Mi c . This figure shows that all 
three components must be substellar. 
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